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Abstract 

We consider the problem of neutral inclusions for two-dimensional conductivity and anti- 
plane elasticity. The neutral inclusion, when inserted in a matrix having a uniform field, does 
not disturb the field outside the inclusion. The inclusion consists of a core and a shell. We 
show that if the inclusion is neutral to two linearly independent fields, then the core and the 
shell are confocal ellipses. 

1 Introduction 
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. The purpose of this paper is to prove that if a coated inclusion in two dimensions is neutral to 

' multiple uniform fields, then the core and the shell are confocal ellipses. 

We begin by defining the neutral inclusions. Let D and fl be bounded domains with Lipschitz 
boundaries in such that D C fl so that D is the core and n \ D is the shell. The conductivity 
is CTc in the core, as in the shell, and am in the matrix (M^ \ ^^)- So the conductivity distribution 
' is given by 

O ■ = ^cX{D) + asx{^ \ D) + araX{^^ \ ^) 

where x is the characteristic function. The conductivities CTc and ag are assumed to be isotropic 
, (scalar), but we allow (t,„ to be anisotropic, i.e., a positive definite symmetric constant matrix. 

For a given function h with V • a„i^h = in M^, we consider 



m 

^ . I V • aVu = in m , 

y—{ • \ u{x) - h{x) = 0{\x\^^) as \x\ ^ oo, ^ ' ^ 

This problem can be regarded as a conductivity problem or an anti-plane elasticity problem. 
I The inclusion ft (or a) is said to be neutral to the field — Vft- if the solution u to (jl.ip satisfies 

i u{x) — h{x) ~ in M? \ Vl. So the neutral inclusion to the field — V/i does not perturb the field 

outside the inclusion. In the imaging point of view, it means that the field — V/i can not probe the 
neutral inclusion. 

A particular interest lies in the inclusions neutral to uniform fields, i.e., h(x) — a ■ x for some 
constant vector a. If D and are concentric disks, say 13 = { < ri } and Vl = {\x\ < r2 }, 
and am is isotropic, then one can easily see that U, is neutral to all uniform fields if the following 
relation holds: 

(cTs + <Jc){<Jm ~ Os) + f[<7s - 0-c)(o-m + a g) = (1.2) 

where f = rl/r^ (the volume fraction) . Much interest in neutral inclusions (to uniform fields) was 
aroused by the work of Hashin |8j. He showed that since insertion of neutral inclusions does not 
perturb the outside uniform field, the effective conductivity of the assemblage filled with coated 
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inclusions of many different scales is am satisfying (jl.2p . It is also proved that this effective 
conductivity is a bound of the Hashin-Shtrikman bounds of the effective conductivity of arbitrary 
two phase composite. We refer to a book of Milton [15' for development on neutral inclusions in 
relation to theory of composites. 

Another interest in neutral inclusions has aroused in relation to imaging and invisibility cloaking 
by transformation optics. In this regard, we first observe that in general the solution u to (jl.ip 
satisfies u{x) — h(x) = 0{\x\~^) as |a;| — > oo. But, if the inclusion is neutral to all uniform 
fields, then the linear part of h is unperturbed and one can show using multi-polar expansions 
that u{x) — h(x) — 0(|x|~^) as \x\ — > oo for any h (not necessarily linear). It means that it is 
harder to probe the neutral inclusions using u — h. Recently, Ammari et al !3_ extend the idea 
of neutral inclusions to construct multi-coated circular structures which are neutral not only to 
uniform fields but also to fields of higher order up to N for a given integer N (which is called the 
GPT- vanishing structure, GPT for generalized polarization tensor), so that the solution u to (|l.ip 
satisfies u{x) — h{x) — 0(|a;|^^"^) as |a;| oo for any h. This structure has a strong connection 
to the cloaking by transformation optics. The transformation optics proposed by Pendry et al 
transforms a punctured disk (sphere) to an annulus to achieve perfect cloaking (the same transform 
was used to show non- uniqueness of the Calderon's problem by Greenleaf et al [7]). Kohn et al 
[Tl] showed that if one transforms a disk with a small hole, then one can avoid singularities of the 
conductivity which occur on the inner boundary of the annulus and achieve near-cloaking instead of 
perfect cloaking. In [3; it is shown that if we coat the hole by the GPT-vanishing structure before 
transformation, then near-cloaking is dramatically enhanced. See also [1] HI [5j [13] for further 
development to Helmholtz and Maxwell's equations. 

There have been some work on neutral inclusions of general shape. Non-elliptic inclusions 
neutral to a single uniform field have been constructed by Milton and Serkov |16) when the con- 
ductivity CTc of the core is either or oo, and by Jarczyk and Mityushev [3] when CTc is finite. In 
|16) . it is also proved that if an inclusion is neutral to all uniform field (or equivalently, to two 
linearly independent uniform fields), then the core and shell are confocal ellipses, when CTc is 
or oo. That the confocal ellipses (and ellipsoids) are neutral to all uniform fields was proved by 
Kerker [H]- 

The purpose of this paper is to prove that confocal ellipses are the only inclusions which are 
neutral to two fields even when Cc is finite. We emphasize that the method of [16] can not be 
applied to the case when CTc is finite. There the conformal mapping from the shell (fl \ D) onto 
an annulus was used. But, there is no conformal mapping from VL onto a disk which maps Z? to a 
concentric disk, except for some very special cases. 

More precisely we prove the following theorem. We emphasize that the theorem holds for any 
conductivity <Tc of the core (finite, 0, or oo). 

Theorem 1.1. Let D and Vl be bounded domains with Lipschitz boundaries in such that D C fl. 
If n is neutral to — Vxj for j — 1,2, then D and il are confocal ellipses. 

The key observation in proving Theorem 1 1.1 1 is that if the inclusion is neutral to two fields, then 
the fields inside the core is also uniform. Using this fact, we are able to set up a free boundary 
value problem. We then use a conformal mapping to show that the solution to the free boundary 
value problem is a pair of confocal ellipses. 

This paper is organized as follows: In the next section we review layer potential representation 
of the solution to (|1.1|) . In section 3, we show that if the inclusion is neutral to two fields, then the 
fields inside the core is uniform, and derive the free boundary value problem. In the last section 
we show that the solution to the free boundary value problem is confocal ellipses, and hence prove 
Theorem O 
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2 Layer potential representations of solutions 



Let B be a bounded domain in with the Lipschitz boundary. The single layer potential Sqb 
of a density function (p e L^{dB) is defined by 



SdBM{x) := — 

ZTT 



dB 



\i\\x - y\ip{y) da{y) , x G 



It satisfies the jump relation 



d 



[lp]{x), x&dB 



where the operator ICgs is defined by 

ICqb[(p]{x) = / 



\x - yY 



(2.1) 



(2.2) 



(2.3) 



and is its L^-adjoint. The subscripts + and — indicate the limits from outside and inside B, 
respectively. 

Suppose that Um is isotropic. It is known (see, for example, [2j) that there is a pair {ip,ip) € 
L§ ((?£)) X L§(9ri) (the subscript indicates the mean value zero) such that the solution u to (jl.ip 
can be represented as 



u{x) ^h{x)+ Sqd Vp\ {x) + Son (x) , x S 



(2.4) 



The transmission conditions (continuity of the potential and the flux) on the interfaces dD and 
dfl are equivalent to 



{xi-ic*oj,M 



d 



dvoD 



d 



dvoQ. 



SodM + - IC*gn)bJ-'] 



dh 

dvoD 
dh 

dvon 



on dD, 
on dVl, 



SoiAi^]^^— ondD, (2.5) 

(2.6) 

^ = 7r/ 7. V- (2.7) 

Z[ac — Crs) l{<7s—CFm) 

Here and throughout this paper Ouoo denotes the outward normal derivative on dD. The system 
of integral equations (|2.5p and (|2.6p has a unique solution (see [2|). We emphasize that the 



where 



representation (|2.4p is valid even if cjc is or cx), and the solution u is defined in Z?. If (Tc = then 
A = —1/2, and if CTc = oo then A = 1/2. 
Note from the jump relation (|2.2p that 
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dvoD 

du 



Sod M 
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■Sol 



+ dvoD 
{h + Sod[v] + SonbiA 
du 



dvoD 



[h + SooVp] +So^ 



duoD 

cr, 



- 1 



dvoD 
du 



dvoD 



So, we have 



2 du 
2A — 1 duoD 



(2.8) 
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Similarly one can show that 



^= 1- — 



du 



du 



2/i + 1 dvQn 



So, if a is neutral to the field — V/i, then u = h inM? \ and hence we have 



dh 



2^ + 1 dvQn ' 

3 Uniformity of the field in the core 

Suppose that a is neutral to both — Vxi and — Va:;2. After a rotation if necessary, we may 
that (Tm is given by 



(2.9) 



(2.10) 



(3.1) 



Let Uj, j = 1, 2, be the solution to (II. ip with h{x) = xj. Then uj is the solution to (|l.ip with 
replaced by cr^j. So, we may represent it as 



Uj{x) = Xj + SdD[(Pj]ix) + Sdfi[ijjj]{x). 
Let n = (711,712) denote the outward unit normal to dQ or dD and let 



2((Ts ~ ai 



Then we have from (j2.5p and (j2.6p 



dvdD 



J = 1,2. 



^aniV'i] = nj on aZ), 



— — Sod Vpj] + - K:*q^)[^j] = Uj on dn, 



and from ((2T0)) 



2^, 



— rij on dVl, j — 1, 2. 



(3.2) 
(3.3) 

(3.4) 
(3.5) 

(3.6) 



m 
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Let Ui be the harmonic conjugate of ui in each connected component of \ [dD U dVl). We 
may choose Ui = X2 in \ f2. We emphasize that even though i7 \ D is not simply connected, 
a harmonic conjugate of ui exists there. In fact, a harmonic conjugate of xi + iSar2[V'i] exists 
fi, and since ipi G L^{dD) and SdoYpilix) — 0(\x\~^), a harmonic conjugate of SqdYpi] exists 
M2 \ U. Define w by 

' al,ui inR2\n, 

(TsUj^ + ci in J7 \ D, (3-7) 
acUi + C2 in D. 

Here constants ci and C2 are chosen so that v is continuous across (917 and dD. Then v satisfies 



where 



f V • (aVv) = 

\v{x) - X2 = 



in M2 \ 

in R2 \ n, 
in 17 \ 
in D. 



(3.8) 



(3.9) 
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So cr is neutral to —Vx2- 
We represent v as 

v{x) = X2 +Sod[^3]{x) + Soni^PaKx). 



Since 



the pair of potential (ipsjips) satisfies 



(-XI-IC*go)[^3]--^Sonm=n2 oiidD, (3.10) 

OVdD 

d 

SaD[(P3] + {-fJ-iI - IC*Qn)[^-'3]^n2 ondfl, (3.11) 



and 



^3 = ^2^-- (3-12) 



It follows from dHU) and ((XT^ that 

2/i2 + 1 

and since Sooi^j] + Sgnl^pj] = outside 17 for j — 2, 3, we have 

S9D[^2]ix) - -5ao[V^2](:^) - -^^^^Sanmix) = ^^^^az^bsK^), e M.^\n. (3.13) 
Since 5aD[<i52] and Sanies] are harmonic in \ I?, we have 



SaD[^2Kx) ^ —^^^—SaD[^3Kx), xeR''\D. 

Ziln + i 



Note that if Saolf] — outside D, then </? = 0. So we conclude 



¥^2 = — r^3- (3.14) 

2/X2 + 1 



We now see that p.lOp can be written as 



-A/ - /Ca£,)[(^2] - T7^'5aa[''/'2] = ,^/ "2 on 91). 
dl/gc 2/X2 + 1 



By comparing this formula with p.4p . we deduce 



We then have from (12.8 



du2 
dv 



2A-1 (2A- 1)(mi +Ai2) 



"^^ ^ 2A(2M2 + 1) 



and hence 

(2A-l)(^i+M2) _^ ^„ ,„^„. 

2A(2/.2 + l) + ^'-''^ 
for some constant C. So — Vm2 is uniform in D. 



By substituting (|3.6p and p.l5|) into (|3.2p . we have 



It then follows from p.l6p that 
(Mi + ^J'2)S^D[n2]{x) + 2XSdn[n2]ix) 



A(2^2 + 1) 



"2 



(a;) + Sdn 



2Ai2 + 1 



"2 



(a;). 



(3.17) 



0, 



- Ai2 - 1) - 



/-il + M2 



X2 + C2, X G D. 

(3.18) 



By the exactly same argument with switched roles of ui and U2, one can show that 

0, 

Ml + M2 



(mi + /^2)'5an[7ii](x) + 2XSon[ni]{x) = 



A(m2 - Ail - 1) - 



Xl + Cl, 



Let No{x) be the Newtonian potential on £>, i.e., 

Nd{x) 7^ / Inlx-yldy, xeR2 

Then, we have 

d 

S9D[nj]{x) = Ar^(x), X e M^. 

It then follows from (|3TT8l) and (|3?T9| that 



X E D. 

(3.19) 

(3.20) 
(3.21) 



(Ail +M2)A^D(a^) +2AiVo(a:) = 
where Ci and C2 are constants and 

di = 
do = 



dixj + d2xl - {cixi + C2X2) + C2 



2X{fii - A^2 + 1) + ipi + M2) 
4 

2A(/^2 - Ail + 1) + (/-ii + M2) 



x€R'^\n, 
X e D, 



We may assume Ci = C2 = by translating D and O. Observe that (/ii + ^2)Nd{x) + 2XNq{x) 
behaves as (mi+M2)PI+2A|o| 

a;[ as |x| — 00, where \D\ denotes the area of D. So, we have 



+ ii2)\D\ + 2X\n\ = 0, 

2A 



-/ 



(3.22) 



Ml + M2 

where / = |Z3|/|ri| (the volume fraction). It is worth mentioning that this relation shows that A 
and Ml + M2 have opposite signs. We finally have 



Nd{x) - fNnix) = <! 1-/(1 + ^1-^2) 2 ^ 1-/(1+M2- Ml) 



-xi + C2 



xeD. 



(3.23) 

It is worth mentioning that if the Newtonian potential of a simply connected domain is quadratic 
inside the domain, then the domain is an ellipse. This fact was proved by Dive |6] and Nikliborc |18j . 
and is the key ingredient in resolution of the Eshelby's conjecture and the Polya-Szego conjecture 
by Kang and Milton [H] (see also [I^). In the next section we show that p.23p implies that D 
and n are confocal ellipses. 
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4 The free boundary value problem 

Let 

w{x) -.^ ^\x\^ + 2{NDix) ~ fNn{x)), xen\D. 

Note that ANfi{x) = 1 for x G 17, ANfi{x) — ioi x ^ fl, and VNq is continuous across dil. So, 
by (|3.23p w satisfies 

' Aw = mn\D, 

Vw = fx on dfl, (4.1) 

Vw ^ X + {fii ~ ^2){xi,-X2) ondD. 



We now show that if the problem ()4.ip admits a solution for some /, then D and Q are confocal 
ellipses. Let z — xi + ix2 and let w be a solution to (|4.ip . Define 

, , „ ^ / X dw . dw 
g{z) ..^2^w{z) = — ~^ — . 

Then g is holomorphic in \ Z?, g{z) = /z on and g{z) = (/ii — ij,2)z + z on dD. Let ^(C) be 
the conformal mapping from A — {( : 1 < \(\ < rg} for some ro onto \ Z). Since 911 and are 
Lipschitz, $ is continuous on A. Let := g(<&(C))- Then /i is holomorphic in A and satisfies 
/J(C) = fW) on Id = ro and /^(C) = ¥(0 + (a^i - M2)1>(C) on [CI = 1. 
Suppose that $ admits the Laurent series expansion 

oo 

^(C)^ E (4.2) 

n— — OO 

Then we have 

' oo 

/ E —.r^^'X" on|C|=ro, 

n— — OO 

oo oo 

E a^C + if^i - ^^2) J2 on|Cl = l. 



HO 



(4.3) 



Tl — — oo 



Since the Laurent series expansion of h is unique, we have 

famro^"" = + (^i - /^2)a„ (4.4) 
for any integer n. Replacing n with — n and taking complex conjugates we also have 

fanrl"" = a„ + (a<i - Ai2)a^- (4.5) 

Above two identities imply 

[(1 - /r„-2")(l - /rg") - (^1 - ^,2r]an = (4.6) 

for any integer n. Since (1 — frQ^"'){l — /rg") takes different values for different positive integers n, 
we know that there are only one positive n with an 0. By (|4.4p . the possibly nonzero coefficients 
of $(w) are a„ and a-„ and the univalence of $ implies n — 1. Finally, dD and 517 are confocal 
ellipses since they are images of concentric circles by the map $(C) — aiC + o,-iC~^- 
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